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Reactions on surfaces play an important role in many tech- 
nological applications. Since these processes are often rather 
complex, one tries to understand single steps of these com- 
plicated reactions by investigating simpler system. In par- 
ticular the hydrogen dissociation on surfaces serves as such 
a model system. There has been much progress in recent 
years in the theoretical description of reactions on surfaces 
by high-dimensional dynamics simulations on potential en- 
ergy surfaces which are derived from ab initio total energy 
calculations. In this brief review I will focus on the hydrogen 
dissociation on the clean and sulfur-covered Pd(100) surface. 
These calculations established the importance of dynamical 
concepts like the steering effect. The electronic structure cal- 
culations allow furthermore the determination of the factors 
that determine the reactivity of a particular surface. This will 
be demonstrated for the poisoning of hydrogen dissociation by 
sulfur adsorption on the Pd(100) surface. In addition, quan- 
tum effects in the dynamics can be assessed by comparing 
the results of classical with quantum dynamical calculations 
on the same potential energy surface. 

68.35. Ja, 82.20.Kh, 82.65.Pa 



I. INTRODUCTION 

Reactions on surfaces play an important role in many 
technological applications like the heterogenous catalysis, 
growth of semiconductor devices, corrosion and lubrica- 
tion of mechanical parts, or hydrogen storage in metals. 
In spite of this significance the understanding of the mi- 
croscopic details of these reactions is still rather incom- 
plete. Of particular importance are processes in which 
chemical bonds of molecules are breaking due to the pres- 
ence of a substrate because there processes represent the 
first elementary step in, e.g., heterogeneous catalysis or 
corrosion. Often this is the rate-limiting step, for ex- 
ample in the ammonia synthesis. Hydrogen dissociation 
on metal surfaces has become the model system for the 
bond-breaking process on surfaces in the last years be- 
cause it can be studied in detail experimentally as well 
as theoretically. In particular in the theoretical descrip- 
tion there has been much progress recently due to the 
improvement of computer power and the development of 
efficient algorithms. It has become possible to map out 
detailed potential energy surfaces of the dissociation of 
hydrogen on metal surfaces by density functional theory 
calculations Jl|-§|. The availability of high-dimensional 
reliable potential energy surfaces has challenged the dy- 



namics community to improve their methods in order 
to perform high-dimensional dynamical studies on these 
potentials. Now quantum studies of the dissociation of 
hydrogen on surfaces are possible in which all six de- 
grees of freedom of the molecule are treated dynamically 
|i] fL3|| , In this brief review I will illustrate this progress 
by focusing on the hydrogen dissociation on the clean 
and sulfur-covered palladium surface. 

Hydrogen is the simplest molecule which makes it ac- 
cessible to a relatively complete theoretical treatment. At 
the same time hydrogen is also well-suited for performing 
experiments which allows a fruitful interaction between 
theory and experiment. I will show that general concepts 
relating to the reactivity of surfaces as well as to dynam- 
ical reaction mechanisms can be deduced from the de- 
tailed comparison of theoretical and experimental results 
of the hydrogen dissociation of metal surfaces. These 
concepts are applicable to any reaction system making 
hydrogen the ideal candidate for studying reactions on 
surfaces. 

II. GENERAL CONCEPTS IN THE 
ADSORPTION DYNAMICS AT SURFACES 

The sticking or adsorption probability is defined as 
the fraction of atoms or molecules impinging on a sur- 
face that are not scattered back, i.e. that remain on the 
surface. It should be noted here that there is no unam- 
biguous definition of the sticking probability because for 
surfaces with non-zero temperature every adsorbed par- 
ticle will sooner or later desorb again. Hence the sticking 
probability depends on the time-scale of the required res- 
idence time on the surface. 

Atomic adsorption is often very efficient. Hydrogen 
atoms, e.g, stick at metal surfaces jl4j and semiconductor 
surfaces |l5| ] with a probability of order unity. However, 
dissociative adsorption probabilities can differ by many 
orders of magnitude. Whereas the sticking probability of 
hydrogen molecules on many transition metal surfaces is 
about 0.5 |l4|,[l6|], at room temperature the dissociation 
probability of H 2 /Si is only 10 -8 @, and for N 2 /Ru it is 
even as low as 10~ 13 |18) . The investigation of processes 
that occur within such a wide range of probabilities rep- 
resents of course a great challenge for the theory as well 
as for the experiment. 

There is one fundamental difference between atomic 
and molecular adsorption on the one side and dissocia- 
tive adsorption on the other side that is very important 
for the theoretical description of these processes. Here 
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FIG. 1. Atomic adsorption probability for Xe and Ar on 
Pt(lll) and molecular adsorption probability of N2/W(100). 
These examples are taken from the textbook by Zangwill J19| . 
The inset illustrates the adsorption process. 

with molecular adsorption a sticking process is meant in 
which the molecule stays intact on the surface. This fun- 
damental difference will be illustrated in the following. 
In atomic and molecular adsorption it is crucial that the 
impinging particles transfer their kinetic energy to the 
surface, otherwise they would be scattered back into the 
gas phase. If -Pe(c) is the probability that an incoming 
particle with kinetic energy E will transfer the energy 
e to the surface, then the atomic or molecular sticking 
probability can be expressed as 

poo 

S(E) = / P E (e)de, (1) 

J E 

i.e., it corresponds to the fraction of particles that trans- 
fer more energy to the surface than their initial kinetic 
energy. This excess energy has to be transferred to sub- 
strate excitation, i.e., either phonons or electron-hole 
pairs. Hence any theoretical description of atomic or 
molecular adsorption has to consider dissipation to the 
continous excitation spectrum of the substrate. In Fig. |l| 
sticking probabilities for atomic and molecular adsorp- 
tion as a function of the initial kinetic energy are shown 
that correspond indeed to textbook examples Jl9| . These 
curves show a typical behavior, namely the decrease of 
the sticking probability with increasing kinetic energy. 
This is due to the fact that the energy transfer to the 
surface becomes less efficient at higher kinetic energies. 
Of course, the higher the kinetic energy is, the more en- 
ergy is transfered to the surface. But the fraction of 
particles that loose more energy than their initial kinetic 
energy becomes smaller at higher kinetic energy. There 
is still more interesting physics in these sticking probabil- 
ities. For example, at low kinetic energies classically the 



FIG. 2. Dissociative adsorption probability versus kinetic 
energy of H2/Cu(lll) for molecules initially in the vibrational 
ground state. Solid line: results of five-dimensional calcula- 
tions in which the molecular axis was kept parallel to the sur- 
face (from ref. |2l|). Dashed line: Experimental curve (from 
ref. P2j). The inset illustrates the dissociation process. 

sticking probability should become unity if there is no 
barrier before the adsorption well. Every impinging par- 
ticle transfers energy to the substrate so that in the limit 
of zero initial kinetic energy all particles will stick. Quan- 
tum mechanically, however, there is a non-zero probabil- 
ity for elastic scattering at the surface so that the sticking 
probabilities become less than unity in the zero-energy 
limit pop . In Fig. [t] these quantum effects at low energies 
are evident in the sticking probability of the light noble 
gas argon on Pt(lll) compared to the sticking probabil- 
ity of the heavier noble gas xenon on the same surface. 

Now in the case of dissociative adsorption there is an- 
other channel for energy transfer, which is the conver- 
sion of the kinetic and internal energy of the molecule 
into translational energy of the atomic fragments on the 
surface relative to each other. This represents the funda- 
mental difference to atomic or molecular adsorption. It is 
true that eventually the atomic fragments will also dissi- 
pate their kinetic energy and come to rest at the surface. 
However, especially in the case of light molecules like hy- 
drogen dissociating on metal surfaces the energy transfer 
to the substrate is very small due to the large mass mis- 
match. Whether a molecule sticks on the surface or not 
is almost entirely determined by the bond-breaking pro- 
cess for which the energy transfer to the substrate can be 
neglected. This makes it possible to describe the dissocia- 
tive adsorption process within low-dimensional potential 
energy surfaces neglecting the surface degrees of freedom 
if furthermore no substantial surface rearrangement upon 
adsorption occurs, as it is usually case in the dissociative 
adsorption on close-packed metal surfaces. Fig. shows 



2 



the dissociative adsorption probability of a system which 
also corresponds to a textbook example, namely the dis- 
sociative adsorption of H2 on Cu(lll). In this system 
the dissociation is hindered by a noticeable barrier so 
that the dependence of the sticking probability on the 
kinetic energy exhibits a behavior typical for activated 
systems @>p|. 



III. 



DISSOCIATIVE ADSORPTION AT A 
TRANSITION METAL SURFACE 



Transition metal surfaces are usually very reactive as 
fas as hydrogen dissociation is concerned |]l6| : p3| -^5| . In 
Fig. H the results of molecular beam experiments of Ren- 
dulic, Anger and Winkler [H| and of Rettner and Auer- 
bach (26) for the dissociative adsorption of H2 on Pd(100) 
are shown. At low kinetic energies these experiments 
yield a sticking probability above 0.5. But even more 
interestingly, the sticking probability initially decreases 
with increasing kinetic energy. This is reminiscent of 
the dependence of the sticking probability on the kinetic 
energy in atomic or molecular adsorption illustrated in 
Fig. |l| Therefore for a long time it was believed that 
such an initially decreasing sticking probability in disso- 
ciative adsorption is a signature of the so-called precursor 
mechanism p7| . In this mechanism the molecule does not 
directly dissociate but is first trapped molecularly in a 
precursor state from which it then dissociates. This trap- 
ping probability decreases with increasing kinetic energy 
and thus determines the sticking at low kinetic energies. 

Wilke and SchefHer have performed density-functional 
theory (DFT) calculations of the interaction of H2 with 
Pd(100) in order to elucidate the dissociation process. In 
Fig. [| two so-called elbow plots are shown. They repre- 
sent a two-dimensional cut through the potential energy 
surface (PES) of Bi2/Pd(100) in which the orientation of 
the molecule, its center-of-mass lateral coordinates and 
the substrate are kept fixed. The molecule is oriented 
parallel to the surface, and the PES is plotted as a func- 
tion of the center-of-mass distance of the molecule from 
the surface Z and the interatomic distance drr-H- Fig. ||a 
demonstrates that the dissociation of H2 on Pd(100) is 
non- activated, i.e., there are reaction paths towards dis- 
sociative adsorption with no energy barrier. The major- 
ity of reaction pathways, however, is hindered by barriers 
p8| . Furthermore, in these calculations no molecular ad- 
sorption state has been found. It looks like there is such 
a well in Fig. |4j). However, the detailed DFT study of 
the PES has shown that this apparent well does not cor- 
respond to a local minimum of the PES, it is rather a 
saddle point in the multi-dimensional PES. 

An analytical representation of this ab initio PES has 
been used for a quantum dynamical study in which all 
six hydrogen degrees of freedom were taken into account 
explicitly while the substrate was kept fixed p8|. The 
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FIG. 3. Sticking probability versus kinetic energy for a hy- 
drogen beam under normal incidence on a Pd(100) surface. 
Theory: six-dimensional results for H2 molecules initially in 
the rotational and vibrational ground state (dashed line) and 
with an initial rotational and energy distribution adequate 
for molecular beam experiments (solid line) || . H2 molecular 
beam adsorption experiment under normal incidence (Ren- 
dulic et al. |l6|): circles; H2 effusive beam scattering exper- 
iment with an incident angle of of 6i = 15° (Rettner and 
Auerbach |Eq|): long-dashed line. 



results are also plotted in Fig. |3|. The sticking curve for a 
monoenergetic beam initially in the vibrational and rota- 
tional ground state shows a strong oscillatory structure 
which will be discussed below. An experimental molecu- 
lar beam, however, does not correspond to a monoener- 
getic beam in one specific quantum state. If one assumes 
an energy spread and a distribution of internal molecu- 
lar states typical for a beam experiment, the oscillations 
are almost entirely smoothed out in the 6D quantum re- 
sults (solid line in fig. |3|). The results corresponding to 
the beam simulation agree with the experimental results 
semi-quantitatively. More importantly, they reproduce 
the general trend found in the experiment, namely the 
initial decrease of the sticking probability as a function 
of the kinetic energy followed by an increase at higher 
energies. Now in the ab initio PES there is no molecu- 
lar adsorption state, furthermore in the 6D quantum dy- 
namical calculation no energy transfer to the substrate 
is considered. Hence the precursor mechanism cannot be 
operative in the simulation. So what is the reason for the 
initial decrease of the sticking probability? 

Since energy transfer cannot play a crucial role in the 
adsorption process, this decrease in the sticking probabil- 
ity has to be caused by a purely dynamical effect, namely 
the steering effect [§,|§j32 



Although the majority of 
pathways to dissociative adsorption has non-vanishing 
barriers with a rather broad distribution of heights and 
positions, slow molecules can be very efficiently steered to 
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FIG. 4. Contour plots of the PES along a two 
two-dimensional cuts through the six-dimensional coordinate 
space of H2/Pd(100), so-called elbow plots, determined by 
GGA calculations . The coordinates in the figure are the 
H2 center-of-mass distance from the surface Z and the H-H 
interatomic distance d. The dissociation process in the Zd 
plane is illustrated in the inset. The lateral H2 center-of-mass 
coordinates in the surface unit cell and the orientation of the 
molecular axis, i.e. the coordinates X, Y, 0, and <j> are kept 
fixed for each 2D cut and depicted above the elbow plots. En- 
ergies are in eV per H2 molecule. The contour spacing in a) 
is 0.1 eV, while it is 0.05 eV in b). 

non- activated pathways towards dissociative adsorption 
by the attractive forces of the potential. This mechanism 
becomes less effective at higher kinetic energies where the 
molecules are too fast to be focused into favourable con- 
figurations towards dissociative adsorption. If the kinetic 
energy is further increased, the molecules will eventually 
have enough energy to directly traverse the barrier region 
leading to the final rise in the sticking probability. 

In order to illustrate the steering effects, we use the re- 
sults of classical molecular dynamics calculations which 
have been performed on exactly the same PES as the 
quantum dynamical calculations J33|. In these classical 
calculations significant differences in the sticking proba- 
bility compared to the quantum results have been found 
which are mainly due to zero-point effects. The steer- 




FIG. 5. Snapshots of classical trajectories of hydrogen 
molecules impinging on a Pd(100) surface. The initial con- 
ditions are chosen in such a way that the trajectories are re- 
stricted to the zz-plane. Left trajectory: initial kinetic energy 
Ei = 0.01 eV. Right trajectory: same initial conditions as in 
the left trajectory except that the molecule has a higher ki- 
netic energy of 0.12 eV. 

ing effect, however, is a general mechanism operative in 
quantum as well as in classical dynamics. I will there- 
fore use snapshots of two typical trajectories in order to 
illustrate the dynamical mechanism responsible for the 
initial decrease of the sticking probaiblity. These trajec- 
tories are plotted in Fig. ||. The initial conditions are 
chosen in such a way that the trajectories are restricted 
to the xz-plane. The left trajectory demonstrates why 
the sticking probability is so large at low kinetic energies 
due to the steering effect. The incident kinetic energy is 
Ei — 0.01 eV. In this particular example the molecular 
axis is initially almost perpendicular to the surface. In 
such a configuration the molecule cannot dissociate at 
the surface. But the molecule is so slow that the forces 
acting upon it can reorient the molecule. It is turned 
parallel to the surface and then follows a non-activated 
path towards dissociative adsorption. This shows how 
molecule with unfavorable initial conditions can still dis- 
sociate due to very efficient steering towards favorable 
configurations. 

This process becomes less effective at higher kinetic 
energies, which is demonstrated with the right trajec- 
tory in Fig. |^. The initial conditions are the same as for 
the left trajectory, except for the higher kinetic energy of 
0.12 eV. Of course the same forces act upon the molecule, 
and due to the anisotropy of the PES the molecule also 
starts to rotate to a configuration parallel to the surface. 
However, now the molecule is too fast to finish this rota- 
tion. It hits the repulsive wall of the PES at the surface 
with its molecular axis tilted by about 45° with respect to 
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FIG. 6. Angular distribution of the in-plane and 
out-of-plane of scattering of H2/Pd(100). The initial kinetic 
energy is Ei = 76 meV, the incident angle is 6i = 32° along 
the (011) direction. The molecules are initially in the rota- 
tional ground state ji = 0. Open circles correspond to rota- 
tionally elastic, filled circles to rotationally inelastic diffrac- 
tion. The radii of the circles are proportional to the logarithm 
of the scattering intensity, x denotes the (Oil) direction, y 
the (011) direction. The specular peak is the largest open 
circle (from Ref. Q). 

the surface normal. At the classical turning point there 
is a very rapid rotation corresponding to a flip-flop mo- 
tion, and then the molecule is scattered back into the 
gas-phase rotationally excited. 

Now I like to come back to the strong oscillatory struc- 
ture of the sticking curve for a monoenergetic beam ini- 
tially in one particular quantum state in Fig. |3|. These os- 
cillations are a consequence of the quantum nature of the 
hydrogen particle, in classical calculations they do not 
appear Jl0| , ^3[ . If a quantum particle is interacting with a 
periodic surface, coherent scattering leads to diffraction. 
Such a calculated diffraction pattern is shown in Fig. ^ 
for hydrogen molecules in the ground state scattering at 
a Pd(100) surface with an kinetic energy of Ei — 76 meV 
under an angle of incidence of 9i = 32°. There are rather 
many diffraction peaks since also rotationally inelastic 
diffraction can occur, i.e., scattering in which the rota- 
tional state of the molecule is changed. Still the num- 
ber of diffraction peaks is finite and increases disconti- 
nously with increasing energy. An analysis of the oscilla- 
tory structure of the sticking probability in Fig. || reveals 
that these oscillations can be related to threshold effects 
associated with the opening of new scattering channels 
|54],[35| . These oscillations have not been observed exper- 
imentally yet, although they have been carefully searched 
for [^6 36 1 . They are very sensitive to the symmetry of 
the initial conditions ]34|,|35|]. Furthermore, since Pd(100) 
is a very reactive surface, a large fraction of the incoming 
hydrogen molecules is not scattered back coherently but 



adsorbs dissociatively. These adatoms then disturb the 
periodicity of the surface and thus suppress, in addition 
to already existing surface imperfections like steps and 
vacancies, the coherence of the scattering events. Hence 
the experimental observation of this oscillations actually 
represents a challenging task. 

The dependence of adsorption and desorption on ki- 
netic energy, molecular rotation and orientatio n p| , ^7[ , 
molecular vibration |2^] , ro- vibrational coupling ]38| , an- 
gle of incidence and the rotationally elastic and in- 
elastic diffraction of H 2 /Pd(100) [Q have been studied 
so far by six-dimensional ab initio dynamics calculations 
on the same PES. The results of these calculations have 
been compared to a number of independent experiments 
|]l6| , j39l -|4l}| , and they are at least in semi-quantitative 
agreement with all of these experiments This shows that 
ab initio dynamics calculations are indeed capable of ad- 
equately describing the hydrogen dissociation on transi- 
tion metal surfaces. 



IV. DISSOCIATIVE ADSORPTION AT A 
SULFUR-COVERED TRANSITION METAL 
SURFACE 

The presence of an adsorbate on a surface can pro- 
foundly change the surface reactivity. A well-known ex- 
ample is the reduction of the activity of the car-exhaust 
catalyst by lead. But also adsorbed sulfur "poisons" this 
catalyst. An understanding of the underlying mecha- 
nisms and their consequences on the reaction rates is 
therefore of decisive importance for, e.g., designing bet- 
ter catalysts. Traditionally an "trial and error" approach 
was used to improve the activity of a catalyst by adding 
some substances. On Pd(100) it is experimentally well- 
known that the presence of sulfur leads to a large reduct- 
ing of hydrogen dissociation probability julQ. While 
at the clean surface the dissociation probability is about 
60% for a kinetic energy of Ei = 0.05 eV, at the sulfur- 
covered surface it drops below 1% at the same energy 
& 

Theoretically this problem had only been addressed by 
a small number of studies. These focused either on the 
adsorbate induced change of the density of states (DOS) 
at the Fermi level 
trostatic field 
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43| (45| or on the adlayer induced elec- 
48). Just recently, the poisoning of 
hydrogen dissociation on Pd(100) by sulfur adsorption 
has been the subject of detailed DFT studies |||9 50 1. 

In Fig. |?] I have collected four elbow plots of the hy- 
drogen dissociation on the (2x2) sulfur covered Pd(100) 
surface |50). In contrast to the clean Pd(100) surface, 
the hydrogen dissociation on the sulfur covered surface 
is no longer non-activated. The minimum barrier, which 
is shown in Fig. ^a, has a height of 0.1 eV and corre- 
sponds to a configuration in which the H2 center of mass 
is located above the fourfold hollow site. This is the site 
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FIG. 7. Cuts through the six-dimensional potential energy surface (PES) of H2 dissociation over (2x2)S/Pd(100) at four 
different sites with the molecular axis parallel to the surface: a) at the fourfold hollow site; b) at the bridge site between two 
Pd atoms; c) on top of a Pd atom; d) on top of a S atom. The energy contours, given in eV per molecule, are displayed as a 
function of the H-H distance, (1h-h, and the height Z of the center-of-mass of H2 above the topmost Pd layer. The geometry 
of each dissociation pathway is indicated in the panel above the contour plots. The large open circles are the sulfur atoms, the 
large filled circles are the palladium atoms. 



which is farthest away from the sulfur atoms in the sur- 
face unit cell. Recall that the most favorable reaction 
path on the clean Pd(100) surface corresponds to the H2 
molecule dissociating at the bridge position between two 
Pd atoms (see Fig. [|a). For this approach geometry the 
dissociation at the sulfur covered surface is now hindered 
by a barrier of height 0.15 eV (Fig. 0b). There is a pe- 
culiar local minimum at the dissociation path for this 
configuration when the molecule is still lA above the 
Pd atoms. There are apparently subtle compensating 
effects between the attractive interaction of H2 with the 
Pd atoms and the repulsion originating from the S atoms. 
Figs. 0a and b show that the dissociation is hindered by 
the formation of energy barriers in the entrance channel 
of the potential energy surface, however, the hydrogen 
dissociation is still exothermic, i.e., the poisoning is not 
due to site-blocking. This result is actually at variance 
with measurements of the hydrogen saturation coverage 
as a function of the sulfur coverage J42[. In these ex- 
periments a linear decrease of the hydrogen saturation 
coverage with increasing sulfur coverages was found. At 
a sulfur coverage of 0s = 0.28, which is close to the one 
used in the calculations, hydrogen adsorption should be 
completely suppressed, i.e., there should be no attractive 
sites for hydrogen adsorption any more. A possible ex- 
planation for this apparent contradiction will be given 
below. 



Closer to the sulfur atoms the PES becomes strongly 
repulsive. This is illustrated in Fig. 0c and d. While 
the dissociation path over the Pd on-top position on the 
clean surface is hindered by a barrier of height 0.15 eV 
B (Fig. |^b), the adsorbed sulfur leads to an increase in 
this barrier height to 1.3 eV (Fig. 0c). Directly above 
the sulfur atoms the barrier towards dissociation even 
increases to values of about 2.5 eV for molecules oriented 
parallel to the surface (Fig. 0d). 

The goal of any theoretical study should be to pro- 
vide a qualitative picture that explains the calculated 
results. There are many different ways of illustrating the 
electronic factors that determine the reactivity of a par- 
ticular system (see, e.g., Refs. [S0SEl|,|l). Current 
studies have emphasized that the reactivity of surfaces 
cannot be solely understood by the electronic density of 
states at the Fermi level In order to understand 

the origins of the formation of the small energy barriers 
at the hollow and bridge site and the large energy barri- 
ers at the top sites, we will therefore compare the whole 
relevant DOS for the H2 molecule in these different ge- 
ometries. For a discussion of the reactivity of the clean 
Pd(100) surface I refer to Ref. j|. Here I focus on the 
changes of the density of states induced by the presence 
of sulfur on the surface (Fig. ||) . 

The information provided by the density of states alone 
is often not sufficient to assess the reactivity of a particu- 
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FIG. 8. Density of states (DOS) for a H 2 molecule situated at (a) (Z,dn-n)= (4.03A, 0.75A) and (b) (Z,dn-n)= (1-61A, 
0.75A) above the fourfold hollow site which corresponds to the configuration depicted in Fig. (?]a, and for a H2 molecule situated 
at (c) (7i,&h-h)= (3.38A, 0.75A) above the sulfur atom which corresponds to the configuration depicted in Fig. ^d. Z and 
du-H denote the H2 center-of-mass distance from the surface and the H-H interatomic distance, respectively. Given is the local 
DOS at the H atoms, the S adatoms, the surface Pd atoms, and the bulk Pd atoms. The energies are given in eV. 



lar system. It is also important to know the character of 
the occupied and unoccupied states. For the dissociation 
the occupation of the bonding a g and the anti-bonding <j* 
H2 molecular levels and of the bonding and anti-bonding 
states with respect to the surface-molecule interaction 
are of particular importance. 

Figure ||a shows the DOS when the H2 molecule is still 
far away from the surface above the fourfold hollow site, 
i.e, in the configuration that corresponds to Fig. [?]a. The 
H-H distance d is 0.75 A and the center of mass of the H2 
molecule is 4.03 A above the topmost Pd layer so that 
there is no interaction between the hydrogen molecule 
and the sulfur covered palladium surface. The large peak 
in the sulfur DOS at -13 eV corresponds to the S 3s 
state. The sulfur p orbitals strongly interact with the 
Pd d states, which is evident from the peak in the sulfur 
DOS at the Pd d band edge (at E-E F = -4.8 eV) and from 
the broad band at higher energies which has substantial 
weight close to the Fermi level. The d band at the surface 
Pd atoms is broadened and shifted down somewhat with 
respect to the clean surface due to the interaction with 
the S atoms jl9| . There is one intense peak in the hydro- 
gen DOS at -4.8 eV which corresponds to the a g state. 
This peak is degenerate with the sulfur related bonding 
state at -4.8 eV, this degeneracy, however, is accidental, 
as will become evident immediately. 

The density of states for the molecule at the minimum 
barrier position of Fig. 0a is shown in Fig. ||b. Now the 
<7 g state has shifted down to -7.1 eV while the sulfur 
state at -4.8 eV remains almost unchanged. This indi- 
cates that there is no direct interaction between hydrogen 
and sulfur. It also proofs that the degeneracy between 



these two states in Fig. ||a is accidental. Furthermore, 
we find a broad distribution of hydrogen states with a 
small, but still significant weight below the Fermi level. 
These are states of mainly H2-surface antibonding char- 
acter (l9|,[30| . A comparison with the hydrogen dissocia- 
tion at the clean Pd surface yields that more H2-surfacc 
antibonding states are populated at the sulfur covered 
surface. This is caused by the sulfur induced downshift 
of the Pd <i-band. These H2-surface antibonding states 
lead to a repulsive interaction and thus to the building up 
of the barriers in the entrance channel of the PES . It 
is therefore an indirect interaction between sulfur and hy- 
drogen that is responsible for the barriers at this site. A 
similar picture explains why for example noble metals are 
so unreactive for hydrogen dissociation: The low-lying 
(i-bands of the noble metals cause a downshift and a sub- 
stantial occupation of the antibonding H2-surface states 
resulting in high barriers for hydrogen dissociation fl5l[ . 

The situation is entirely different if the molecule ap- 
proaches the surface above the sulfur atom. This is 
demonstrated in Fig. ||c. The center of mass of the 
H2 molecule is still 3.38 A above the topmost Pd layer, 
but already at this distance the hydrogen and the sulfur 
states strongly couple. The intense peak of the DOS at 
-4.8 eV has split into a sharp bonding state at -6.6 eV 
and a narrow anti-bonding state at -4.0 eV. Thus it is 
a direct interaction of the hydrogen with the sulfur re- 
lated states that causes the high barriers towards hydro- 
gen dissociation close to the sulfur atoms. In conclusion, 
the poisoning of hydrogen dissociation on Pd(100) by ad- 
sorbed sulfur is due to a combination of a indirect effect, 
namely the sulfur-related downshift of the Pd <i-bands 
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resulting in a larger occupation of H2-surface antibond- 
ing states, with a direct repulsive interaction between H2 
and S close to the sulfur atoms. 

In order to assess the dynamical consequences of 
the sulfur adsorption on the hydrogen dissociation, six- 
dimensional dynamical calculations on the analytical rep- 
resentation of the ab initio PES of H 2 at S(2x2)/Pd(100) 
have been performed ]55| . The results of these quantum 
and classical calculations for the H2 dissociative adsorp- 
tion probability as a function of the incident energy are 
compared with the experiment p6| in Fig. ^|. In addition, 
also the integrated barrier distribution Pb(E), 



Pb(E) 



1 



Q{E-E h {9,cj), X, Y)) 



2ttA 

x cos 9d9 d(j) dX dY 



(2) 



is plotted. Here 9 and (f> are the polar and azimuthal ori- 
entation of the molecule, X and Y are the lateral coordi- 
nates of the hydrogen center-of-mass. A is the area of the 
surface unit cell. Each quadruple defines a cut through 
the six-dimensional space (see Fig. for examples), and 
i?b is the minimum energy barrier along such a cut. The 
function 8 is the Heavyside step function. The quantity 
Pb(E) is the fraction of the configuration space for which 
the barrier towards dissociation is less than E. If there 
were no steering effects, Pb(E) would give the classical 
sticking probability, i.e., it corresponds to the sticking 
probability in the classical sudden approximation or the 
so-called "hole model" |§. 

First of all it is evident that the calculated sticking 
probabilities are significantly larger than the experimen- 
tal results. Only the onset of dissociative adsorption at 
E[ « 0.12 eV is reproduced by the calculations. This on- 
set is indeed also in agreement with the experimentally 
measured mean kinetic energy of hydrogen molecules des- 
orbing from sulfur covered Pd(100) [57| ], The question 
arises where these large differences between theory and 
experiment come from. It might be that uncertainties in 
the experimental determination are responsible for the 
difference. The exact sulfur coverage in the experiment 
was not very well characterized. The sulfur adlayer was 
obtained by simply heating up the sample which leads 
to segregation of bulk sulfur at the surface. The sulfur 
coverage was then monitored through the ratio of the 
Auger peaks Si32/Pd33o [^6[ . Since the hydrogen stick- 
ing probability depends sensitively on the sulfur cover- 
age |uj,^g, a small uncertainty in the sulfur coverage 
can have a decisive influence. However, as noted above, 
while the DFT calculations yield that the poisoning is 
caused by the building up of barriers hindering the dis- 
sociation, the vanishing hydrogen saturation coverage for 
roughly a quarter monolayer of adsorbed sulfur 42 sug- 
gests that any attractive adsorption sites for hydrogen 
have disappeared due to the presence of sulfur. These 
seemingly contradicting results and also the discrepancy 
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FIG. 9. Sticking probability versus kinetic energy for a H2 
beam under normal incidence on a S(2x2)/Pd(100) surface. 
Full dots: experiment (from ref. Dashed-dotted line: In- 

tegrated barrier distribution, which corresponds to the stick- 
ing probability in the hole model |56| ; Solid line: Quantum 
mechanical results for molecules initially in the rotational and 
vibrational ground-state; Dashed line: Classical results for 
initially non-rotating and non- vibrating molecules. The inset 
shows the quantum and classical results at low energies. 



between calculated and measured molecular beam stick- 
ing probabilities could be reconciled if subsurface sulfur 
plays an important role for the hydrogen adsorption en- 
ergies. Subsurface sulfur is not considered in the calcula- 
tions but might well be present in the experimental sam- 
ples. The possible influence of subsurface species on re- 
actions at surfaces certainly represents a very interesting 
and important research subject for future investigations. 

Except for this open question, there are further in- 
teresting results obtained by the dynamical calculations. 
The calculated sticking probabilities are not only much 
larger than the experimental ones, they are also much 
larger than what one would expect from the hole model. 
This demonstrates that steering is not only operative 
for potential energy surfaces with non-activated reaction 
paths like for H2/Pd(100), but also for activated systems 
as H 2 /S(2x2)/Pd(100). As Fig. demonstrates, the sul- 
fur covered Pd surface represents a strongly corrugated 
system with barrier heights varying by more than 2 eV 
for molecules with their axis parallel to the surface. And 
the large barriers above the sulfur atoms extend rather 
far into the gas phase (see Fig. |7|d). Molecules with un- 
favorable initial conditions are very effectively reoriented 
to low-barrier sites p5[ . This leads to an enhancement 
of the sticking probability with respect to the hole model 
by a factor of three to four. 

Figure H shows furthermore that the classical molecular 
dynamics calculations over-estimate the sticking proba- 
bility of H 2 at S(2x2)/Pd(100) compared to the quantum 
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results. At small energies below the minimum barrier 
height the quantum calculations still show some disso- 
ciation due to tunneling, as the inset of Fig. ^| reveals, 
whereas the classical results are of course zero. But for 
higher energies the classical sticking probability is up to 
almost 50% larger than the quantum sticking probabili- 
ties. This suppression is also caused by the large corru- 
gation and the anisotropy of the PES. The wave function 
describing the molecule has to pass narrow valleys in the 
PES in the angular and lateral degrees of freedom in or- 
der to dissociate. This leads to a localization of the wave 
function and thereby to the building up of zero-point en- 
ergies which act as additional effective barriers. While 
the vibrational H-H mode becomes softer upon disso- 
ciation so that the zero-point energy in this particular 
mode decreases, for the system H 2 /S(2x2)/Pd(100) this 
decrease is over-compensated by the increase in the zero- 
point energies of the four other modes perpendicular to 
the reaction path, i.e., the sum of all zero-point energies 
increases upon adsorption J55|. Therefore the quantum 
particles experience an effectively higher barrier region 
causing the suppressed sticking probability compared to 
the classical particles. Interestingly enough, if the sum 
of all zero-point energies remains approximately constant 
along the reaction path as in the system H2/Pd(100), 
then these quantum effects almost cancel out [ ]l(i| , |33|] . 

V. CONCLUSIONS 

In this review ab initio studies of reactions on surfaces 
have been presented. In the last years the interaction 
between electronic structure calculations on the one side 
and dynamical calculations on the other side has been 
very fruitful. The availability of high-dimensional re- 
liable potential energy surfaces has challenged the dy- 
namics community to improve their methods in order to 
perform high-dimensional dynamical studies on these po- 
tentials. Now quantum studies of the dissociation of hy- 
drogen on surfaces are possible in which all six degrees of 
freedom of the molecule are treated dynamically. In this 
review I have tried to show that this achievement rep- 
resents an important step forward in our understanding 
of the interaction of molecules with surfaces. Not only 
the quantitative agreement with experiment is improved, 
but also important qualitative concepts emerge from the 
electronic structure calculations as well as from the high- 
dimensional dynamical simulations. These concepts are 
applicable to any reaction system. This represents the 
importance of hydrogen as a model system for studying 
reactions on surfaces. However, the ab initio treatment 
of reactions on surfaces has now matured enough so that 
in the future there will be also more studies on other 
reaction systems like for example the important class of 
oxidation reactions on surfaces. This will be the next 



step towards a full microscopic description of catalytic 
reactions, one of the ultimate goals of surface science. 
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